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FLIGHT SIMULATION OF THE 
LONGITUDINAL MOTIONS 
OF ANOTHER AIRCRAFT 

* SUMMARY 

Simulation of the control -fixed longitudinal dynamic character- 
istics of the A^D-2, a Douglas jet attack aircraft, was attempted 
with use of a variable-stability airplane, Princeton University's 
North American Navion, and with telemetry for data presentation. 

The variable -stability installation was built as a joint project 
for both longitudinal and lateral -directional studies, and although 
the attempt at simulation of the AhD-2 was not successful, it is 
hoped that suitable guide lines have been established for similar 
studies. As a secondary objective, the effects of feeding back 
the airplane responses u,c*, 0 and 0 were investigated. 

The project commenced in October 1959 at Forrestal Research 
Center, Princeton University, Princeton, New Jersey, and was 
completed in May i960. 
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FLIGHT SIMULATION OF THE 
LONGITUDINAL MOTIONS 
OF ANOTHER AIRCRAFT 

INTRODUCTION 

Variable -stability aircraft and their use in flight research have 
received increased emphasis with the advent of sonic and supersonic air- 
craft. To date, much work has been done along these lines, primarily by 
Cornell Aeronautical Laboratory, Inc., and the NACA (now NASA). For 
longitudinal studies in particular, several aircraft including a C-4-5, 
an F-94, a B-2 6 , and a T-33> have been variously configured by Cornell 
Aeronautical Laboratory with electro -hydraulic systems used to modify 
elevator and other auxiliary control surface responses, in order that 
studies could be made of : 

1) aircraft natural modes that are essentially aperiodic and 
convergent; 

2) optimum handling characteristics; 

3) proposed handling characteristics of an aircraft in design stages 

4) inertia effects; 

5 ) phugoid damping . 

Similarly, the NACA has devoted much effort towards such studies, also by 
use of variously configured variable -stability aircraft. Refs. 1 through 
8 apply. ■ 



In this report, simulation of the longitudinal response character- 
istics of the A4D-2, the Douglas attack jet aircraft, was attempted 
with a North American Navion configured as a variable-stability air- 
craft. Only specific response characteristics of the phugoid and short 
period were studied, not pilot opinion or handling qualities, since 
stick-force gradients and the like were not provided. 

The variable -stability installation built in Princeton University's 
North American Navion was constructed around a 3 -axes, stick-steering, 
attitude-and-rate-sensing auto-pilot. Electrical signals produced by 
various transducers were summed in the auto-pilot summing circuits, and 
the control surfaces were moved by electrical servos through the ship's 
control cable system in proportion to these signals. The variable- 
stability installation was used to alter the normal stability derivatives 
of the Navion and to introduce additional or psuedo-stability derivatives, 
thereby changing the aircraft dynamic response characteristics . 

The auto-pilot installation was a joint project in which both longi- 

> 

tudinal and lateral-directional characteristics were studied. The other 
half of the joint project, lateral -directional, is described in Ref. 9* 

All studies were performed at the Forrestal Research Center, Princeton 
University, during the academic year 1959-1960. 
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THEORY AND ANALYSIS 



The normal procedure in the determination of aircraft dynamic motions 
is to solve the simultaneous, constant -coefficient, linear differential 
equations of motion that characterize the transient response of the air- 
craft to a disturbance from some equilibrium flight condition. This 
procedure is outlined in Ref. 10. In the determination of the design 
requirements of the variable-stability installation, the longitudinal 
equations of motion of both the Navion and the A^D-2 were compared as 
outlined below. 



EQUATIONS OF MOTION 

The longitudinal (stick-fixed) equations of motion may be written 
as follows, in non-dimensional time with the angles in radians: 
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Nov let the forcing function, £ e , take the following non-dimensional 
form with angles in radians; 

<S e = K s & s + K l<* d9 + k 3 u + K 4 e (M 

where K^, K^, and hereafter will he referred to as gain 

constants. Equation (4) may also be written as: 
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where £ s , u, <X ,9 and 9 are understood to be perturbation 
quantities . 

Substituting £ into the moment equation : 
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To write the moment equation in real time, p =— - d - Transposing 
and collecting terms : 
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Expressing all the equations in a determinant form in real time: 
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GAIN CONSTANTS 

Solutions to these equations of motion were obtained by assuming 
the solution of the form: 



u = u^e 



K t/r 



C< = <X^e 



x t/r 



9 - Sje * t/r 



Upon substitution of these variables into the equations of motion, 
together with their first and second derivatives, three homogeneous 
algebraic equations resulted. Further, expansion of the resultant 
determinant gave the standard quartic equation in K , implicitly in 
terms of the gain constants, as shown: 
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In terms of modified transfer functions, the numerator terms were 
found to be: 
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Evaluation of the coefficients of the characteristic equation, in 
terms of the Navion stability derivatives as shown in Appendix A, yielded 



the following relations explicit in the gain constants: 

C k = 1 (l 6 ) 
C 3 = 5-30 + 15-75 Kg (17) 
C 2 = 9-38 + 15-74 K x + 32.27 Kg + 15.73 K h (18) 
Ci = .408 + .468 K x + 1.284 Kg - 1.942 K 3 + 32.25 K4 (19) 
C Q = -355 + 1.051 % - 5.81 K 3 + 1.285 K k (20) 



Also as shown in Appendix A, like constants of the characteristic 
equations of the Navion and the A4D-2 were equated in the following 
fashion, and solutions for the gain constants were found: 

( c 424d-2 = ^ c 4^mod. NAV. 

( c 324d_2 = ^ C 3^mod.nav. 

and so on, where the notation "Mod. Nav." refers to Modified Navion. The 
gain constants were found to be as follows : 
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SIGNIFICANCE OF GAIN CONSTANTS 



As may be seen in Equation (6), additional airplane moments were 
added to the longitudinal moment equation. For instance, the C^G^u 
term is implicitly composed of the two terms C^ + ^m&£ K^, where 

c m 6 e K 3 is a moment increment added to the original moment C^. Then 
the individual moments are as follows: 
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Examination of these equations shows the following: 

1) for the Navion, C m = 0, so AC m amounts to a pseudo static 

U1 u “u 

stability derivative; 

2) the same applies to AC^; 

3) since aircraft responses to o< and 0 already exist, A C m ^ 

and AC^q amount to altered stability derivatives and changed 
response characteristics of the aircraft. 

Signs of the gain constants indicate that through supply de- 
stabilizing moments, while supplies a stabilizing moment. A study of 
the magnitudes of the gain constants shows the , following: 

l) K^: for a positive 15 ° increase in angle of attack, the variable- 

,.0 

stability installation must apply .6 of up elevator; 
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2) K^: for a positive pitch rate of 5°/sec., 1.2° up elevator; 

3) K^: for an increase in speed of 12 mph, .3° down elevator; and 

4) k^: f° r a positive 15° increase in pitch angle, .05° down elevator. 
As was shown by Equations (13) through (15), the forcing function 

c m £ e K s & s (which is effectively C m ^ Q £> e ) appears in the numerators, 
but the gain .constants do not so appear. This is clearly the result of 
making changes only to the moment equation, and none to the lift and drag 
equations . The above indicated that with provisions only for modifying 
moments, no simulation could be made of: 

1) amplitude of oscillations in u, oC , 0, and 0 during transient 
responses, and 

2) steady-state values of the above for say an elevator step input 
forcing function. 

It was seen, however, that this would not limit simulation of the transient 
response characteristics damping and frequency. 
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PROCEDURE AND EQUIPMENT 



ANALYTICAL STUDIES 

It was decided that a complete analytical study should precede the 
flight simulation, with particular emphasis on the magnitude of acceleration 
normal to the aircraft's flight path, and some emphasis on the effects of 
the impossibility of simulating aircraft transfer functions of the A4D-2. 

In order to accomplish the above, a computer program utilizing the Goodyear 
Geda Analog Computer was organized. The dynamic response characteristics 
of the A4D-2 were studied first, and damping and frequency of each of the 
two longitudinal response modes, short period and phugoid, were determined 
from the computer and compared with those obtained by analytical methods. 

For comparison, response characteristics of the unmodified Navion were 
also determined, as outlined in Appendix B. Then by use of the computer 
the simulation of the A4D-2 was accomplished, and the effects of vari- 
ations of gain constants were studied. Normal accelerations, primarily 
during the short period, were calculated, also as outlined in Appendix B. 

In addition to the above, studies were made of the effects of the 
gain constants on aircraft response as each individual gain constant was 
fed back with increasing magnitude from a zero value. 

DESIGN OF THE AUTO -PILOT 

Simultaneously with the initial theoretical studies of the variable - 
stability installation, physical design considerations of the auto-pilot 
variable -stability combination were examined. The first step was a flight 
test program to determine control cable forces encountered in the Navion 
during maneuvering flight. Flight test instrumentation used in these 
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flight tests was kept simple to avoid an excessive delay in the construction 
of the auto-pilot. Test instruments included a hand-held spring force gage, 
a glass tube -and -spring accelerometer, and elevator, aileron, and rudder 
control position indicators, which consisted of markings on bulkheads of 
the aircraft close to the respective control and pointers attached to the 
cockpit controls. Since auto-pilot construction had already begun on 
Navion N 91566 , Princeton University's second Navion, N 5 H 3 K, which is 
essentially an identical airplane, was used for the flight tests. 

Flight tests consisted of one flight where measurements were taken of 
control stick force per g, stick force per velocity, and rudder and aileron 
control force per deflection of each. Control stick force per g was measured 
in steady turns, stick force per velocity was measured over a range of 
velocity below and above the trim speed, 122 m.p.h., in 10 m.p.h. increments, 
and rudder and aileron control force per deflection were measured in steady- 
sideslips. One test altitude was used, 3000 ft., and one configuration: 
landing gear and flaps up, canopy closed, cowl flaps closed, mixture rich, 
1900 rpm, and 22 in. of manifold pressure. All tests were performed as 
outlined in Refs. 13 and l4. 

With the Navion on the ground, the hand-held spring force gage was 
used to determine control cable forces encountered during the above flight 
test work. Although measurements of control cable forces were made diffi- 
cult by a cross -control coupling between the aileron and rudder cables, it 
was believed that these forces could be subtracted out with sufficient 
accuracy for design purposes. It was subsequently accomplished in this 
manner. 
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Cable force determination was the linking step to the determination 
of required servo pulley gearing ratios. The two basic considerations in- 
volved in the gearing ratios were whether the given servo-pulley combination, 
with l60 in. lb. design torque, could provide the required cable forces, and 
whether the speed of response of the servos was such that the airplane -auto- 
pilot combination would be dynamically stable. The latter consideration is 
outlined in detail in Ref. 9, and the result of both considerations led to 
the design and construction of a larger servo pulley for the longitudinal 
axis for better speed of response. Existing lateral and directional servo 
pulleys were found to be satisfactory. 

With the help of detailed blueprints and some parts of a Lear, Inc. 
electric auto-pilot specifically designed for the Navion, location of 
separate variable -stability auto-pilot components was decided. With limited 
space available in the Navion, and with the existing arrangement of control 
cables, the logical location of the three -axis servo installation was the 
center of the bulkhead at the after end of the baggage compartment. De- 
sign and construction of the servo support mountings and brackets commenced 

early in November. Design and construction of the electrical and electronic 

% 

components got underway about the same time, as outlined in Ref. 9* 

Angle of attack and sideslip angle sensors were procured from the Naval 
Air Test Center, Patuxent River. By use of a probe, previously designed for 
similar measurements, the angle of attack and sideslip angle’ sensors were 
mounted in tandem, approximately one and one -half chord lengths ahead of 
the leading edge of the left wing, at the tip. The dynamic and static 
pressure probe, also previously designed, was mounted on the right wing-tip, 
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also approximately one and one-half chord lengths ahead of the leading edge; 
the transducer used for pressure (airspeed) measurements was located just 
inside the right wing tip bow. 

An electrical stick and rudder was designed and installed as described 
in Ref. 9 j along with the feed-back gain potentiometers and their mountings. 
The balance of electrical and electronic equipment including telemetry was 
designed and was installed in the baggage compartment; see Ref. 9* 

When the servo and pulley installation had been completed, elevator 
servo cables were fed aft from the servo installation and tied into the ship's 
elevator cables at a small angle, with use of cable tie-ins provided by Lear; 
aileron servo cables were fed forward through a pulley system to the then 
existing aileron follow-up cables, and rudder servo cables were fed aft 
through a single pulley to the ship's rudder cables. Tension on the ship's 
cables was set at 30 lbs., as per Navion maintenance specifications, and 
tension in servo cables was set at 25 lbs . 

Servo clutch slip levels were set such that the safety pilot, with use 
of the ship's controls, could, for example, over-ride the elevator servo of 
the auto-pilot with 1 6 lbs. of control force. In terms of torque limits 
for over-ride, slip levels were set as follows: elevator servo, 160 in. lbs.; 
aileron servo, 100 in. lbs.; rudder servo, 60 in. lbs. These levels were 
designed to permit control of the airplane in the cases of either inadvertent 
engaging of the auto -pilot, or failure of the auto-pilot disengaging system. 

Electrical and mechanical installations were mostly completed by the 
middle of February, and a complete static (ground) test of the servo system 
was inaugurated. The static check included: 
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1) setting the servo limit switches; 

2) with full control deflection, checking the slack sides of 
the control cables for tension; 

3) measuring the control force required to overpower the 
auto -pilot; 

4 ) observing the structural effects on the servo and pulley 
installation, and surrounding mountings, during maximum 
force tests; 

5) observing all cables for excessive cable stretch during 

rapid, large deflections of all controls; 

6) checking the complete aircraft control system for proper 
and safe operation; 

7) inspecting the aircraft structurally after completion of 
the above; and 

8) testing for total weight and aircraft center of gravity 
location. 

Upon completion of the static tests and electrical trouble-shooting, 
at the end of March, the aircraft was test-flown to determine possible 
discrepancies of the aircraft and the auto-pilot installation. After the 
test flights, the Navion was given a somewhat less strenuous but thorough 
static check. 

CALIBRATIONS 

Calibrations consisted essentially of determining first the number of 
degrees of control surface deflection per volt input to the servo, for a 
given servo feed-back gain, and second the magnitude of input signal from 
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each transducer for particular potentiometer settings, e.g., the number 
of volts input from the angle of attack vane at 5° angle of attack, with 
a potentiometer setting of 0.1. The four sensors of the longitudinal 
channel, angle of attack vane, pressure (airspeed) transducer, pitch gyro 
and pitch -rate gyro, were calibrated in this fashion, as was the electrical 
stick. Aircraft response quantities which were to be telemetered, the four 
above plus elevator servo position (elevator position) were calibrated also; 
this amounted to measuring the telemeter d.c. output voltage for an a.c. 
transducer signal input. 

FLIGHT TEST PROGRAM 

The flight test program consisted of three flights which were necessary 
to trouble-shoot communications, telemetry, and auto-pilot electronics, and 
one flight during which test data was taken. Time did not permit further 
testing. 

The configuration tested was as follows: landing gear and flaps up, 

j 

canopy closed, oil cooler doors open, 1850 r.p.tn., and approximately 20 in. 
of manifold pressure. Test altitude was 6500 ft., and airspeed was 110 m.p.h. 
indicated (120 m.p.h. true). The forcing functions utilized in the tests 
consisted of approximately .05 rad. elevator steps or impulses; these were 
provided by the auto-pilot with deflection of the electric stick, which 
had been calibrated to produce approximately .05 rad. of elevator deflection 
for full aft deflection of the stick. Also, during phugoid tests, the wings 
were held level by use of the auto-»pilot rudder controls . 

First, dynamic characteristics of the basic Navion were tested with 
the auto-pilot and the elevator step forcing function; with telemetry, the 
dynamic response characteristics were recorded on magnetic tape by an 
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Ampex 309 -C tape recorder and on time history traces by a Sanborn 150 
recorder. Second, with an arbitrary gain potentiometer setting of .15 
for the u feed-back, and with the auto -pilot and a .05 rad. impulse 
forcing function, the short period and phugoid motions of the modified 
Navion were tested and recorded. In similar fashion, the and 0 

feed-backs weftre tested. After the completion of each test run, the auto- 
pilot was disengaged, the aircraft was flown back to test altitude and 
test airspeed, and the auto-pilot was re-engaged. 

For reasons outlined in later sections, the simulation of the AUD-2 
was not attempted. 
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RESULTS AND DISCUSSION 



GENERAL RESULTS 

The primary objective, that of in-flight simulation of the A4D-2 
with the Navion and its variable -stability installation, was not 
accomplished because of the effects outlined below. The secondary 
objective, which developed during the theoretical analysis of the simu- 
lation, that of an investigation of the effects of variation of the 
feed-back quantities on airplane dynamic response, was accomplished 
through a theoretical analysis, but was not verified by flight test data. 
However, it is felt that the information contained herein concerning the 
variable -stability installation and the attempted simulation should be of 
help in establishing guide lines for similar studies. 

ANALYTICAL RESULTS 

As outlined in some detail in Appendix B, a comparison of the dynamic 
characteristics of the two aircraft, the A4D-2 and the Navion, was made. 
For ease of comparison, these characteristics are repeated here; factor- 
ization of the two characteristic equations gave, in seconds: 
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From the above, it may be seen that the variable-stability installation 
had to be designed to change the dynamic characteristics of the Navion 
in the following way: decrease the frequency of the phugoid and the short 

period modes, decrease the damping of the short period mode, and increase 
the damping of the phugoid mode. As outlined in the section on theoretical 
considerations', such changes were provided for by the use of the feed-backs 
u,c* ,9, and 9, and the related gain constants were found to be of particular 
magnitudes and algebraic signs. 

At this point it should be re-emphasized that the characteristic 
equation of the Navion was made exactly similar to that Df the A4D-2 by 
use of the feed-backs, but, in terms of transfer functions, the numerator 
terms were in no way changed. The transfer functions of the modified 
Navion then consisted of a constant comprised of Navion stability deri- 
vatives, a fourth order denominator term peculiar to the AUD-2, and a first 
or second order numerator term peculiar to the Navion. 

As outlined in Appendix B, the effects of this kind of matching were 
investigated. In a comparison of steady-state values, i.e., the magnitudes 
of the perturbation quantities u,c<,9, and 0 as t-*-a>, the same amount of 
elevator deflection used as a step forcing function produced markedly 
different results. For example, a .05 rad. elevator step produced 
u ss = ,2b6 in the Ato-2, and u ss = 6.15 in the modified Navion. As noted 
in the Appendix, the latter is well beyond the small perturbation theory 
employed in this investigation. Examination of the individual stability 
derivatives involved in the steady state transfer functions showed why the 
steady-state values differed so markedly between the two aircraft. Numer- 
ator terms causing the difference include C m and , and a denominator 
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of the 



term, C^, further complicated the picture; For example, C m 

Navion is approximately four and one-half times the value of that of the 

A4D-2; C t of the Navion is about twice that of the A4D-2; and in the 
L <* 

denominator, C T of the Navion is about half that of the A4D-2. When 

these parameters are multiplied together, as is done in the transfer 

function, the steady-state value of u of the Navion, considering an 

equal amount of elevator step input, should be approximately 20 times 

that of the A4D-2. It should be emphasized that the principal offender 

in the steady-state transfer functions is the derivative C m . In gener- 

oe 

al, the same comments apply to the other perturbation quantities, o< ,9 
and 9, to varying degree. The above results were verified by computer 
studies and are tabulated in Table I. 

As mentioned in EQUIPMENT AND PROCEDURE, computer studies were also 
made of the dynamic characteristics of the A4D-2, the Navion, and the 
modified Navion. Dynamic characteristics such as frequency and damping 
of the oscillatory modes peculiar to the "three" airplanes are also listed 
in Table I, and generally good agreement may be seen, except in the case 
of frequency and damping of the short periods; accuracy in the determination 
of short period characteristics from a time history trace is difficult at 
best. Fig. B1 shows the computer wiring diagram used in the computer simu- 
lation of the A4D-2, and it may be noted that it is a true representation 
of the type of feed-back system used in the Navion variable -stability 
installation . 

Results of these computer studies also included transient values of 
the dynamic responses u,c* ,b and 9, for both the A4D-2 and the modified 
Navion. These results are also tabulated in Table 1. It may be noted 
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that the initial amplitudes of the response quantities to an elevator step 
forcing function of magnitude .05 rad. appear reasonable for the A4D-2, in- 
cluding normal acceleration(g's), but such quantities for the modified 
Navion are not only markedly different, but are well out of the realm of 
the small perturbation theory, as was true in the case of steady-state 
values. For the above reason, emphasis was shifted to the use of an im- 
pulse forcing function; as may be seen in Table I, however, even a .05 rad. 
impulse forcing function produced large initial amplitudes in oscillations 
of the response quantities, including 3*2 to 4.0 g's during the short 
period mode. 

In order to determine the cause of such large differences between the 
two aircraft in amplitudes of oscillations of the response quantities, the 
airplane transfer functions were again examined, with emphasis this time 



on transient response. This procedure is also outlined in Appendix B. Ex- 

u 

Te 



u , . 

ami, nation of the transfer function ^-(S) showed that the gain of the trans- 



fer function is composed of, in the numerator, C m and Cp , and in the 

8 fi a 

3 

denominator, a constant 2, h, and 7" . Evaluation of the gain showed that 
C m and T were the main causes of the large difference in gain between 
the A4D-2 and the modified Navion, viz., the gain of the Navion was about 
four times that of the A4D-2. Examination of the transfer function 

<=e 

showed that the gain of the modified Navion was approximately seven times 
that of the A4D-2, and to further complicate the picture, the damping ratio 
of the second order lead term of the numerator of the Navion transfer 
function is approximately one-third that of the A4D-2, and well below criti- 
cal damping. These factors appreciably altered amplitudes of response. 
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Then it may be seen that C m again altered desired response character- 
istics of the modified Navion, and additionally, the terms h and f took part 
in this alteration. Examination of these terms and others will now follow, 
and later sections will touch on the same subject. It may be noted that a 
reduction of the magnitude of the elevator step or impulse will reduce the 
amplitude of bscillations of the response quantities of the modified Navion, 
and such reduction is linear with decrease in size of the step or impulse. 
However, there is some lowqr limit to this type of reduction; from the pre- 
vious section it may be seen that the .05 rad. elevator step, for example, 
should be reduced by a factor of about 20 to . l4°. Examination of this from 
a practical standpoint, considering cable slop and resolution of the average 
auto-pilot, shows that this is the right direction but the wrong path. 

Examination of the effects of C m indicates that if C mc could be re- 

oq, t be 

duced, response characteristics would be improved; a solution in this direction 
would entail a small, auxiliary surface to be used as a moment -changer vice 
the Navion elevator. This will be discussed in more detail in a later section. 
Reduction of the amplitudes of the oscillations could be effected by increas- 
ing the magnitudes of h and y ; this would indicate that the airplane to be 
simulated should have approximately the same mass and the same moment of 
inertia characteristics about the Y-axis as does the simulator airplane. It 
may be shown, however, that with a 4- feed-back, moment of inertia character- 
istics may also be simulated. Then, in the general case, y should be made 
approximately equal to that of the simulated aircraft. 

In the last direction indicated above, that of using a light air- 
plane to simulate a light airplane, resolution of the auto-pilot again 
becomes part of the picture. This phase was also investigated, and as 
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was expected, in the case that the aircraft and simulator aircraft are 
much alike in terms of stability derivatives and mass characteristics, 
the gain constants, tend to become extremely small; then the 

required magnitudes of the gain constants would tax an auto -pilot which 
had resolution approaching infinity — or characteristics of a pure 
analog computer. 

Because of the above, and other effects outlined in following 
sections, and the very limited time remaining for further investigation, 
emphasis was shifted to the secondary objective mentioned previously: 
a qualitative investigation of the effects of variation of the feed-back 
quantities. This was accomplished by means of root locii plots and was 
verified ty computer studies. Only feed-backs of the algebraic sign de- 
signed for the A^D-2 were studied. Details of the root locii plots may 
be seen in Appendix B, and the general effects of the feed-back quantities 
are summarized in Table II. General comments concerning these effects are 
as follows: the principal effect of feeding -back angle of attack is to 

cause the frequency of the short period mode to decrease and the mode 
tends to become a convergent, aperiodic motion; there is little effect on 
damping of the short period, and frequency and damping of the phugoid. 

Feeding back pitch -rate has the effect of altering dynamic character- 
istics of both the short period and phugoid modes; at very high gain, the 
frequency and damping of the short period mode decrease, and the mode 
tends to become aperiodic, convergent. Further, at low gain the damping 
and frequency of the phugoid increase, at some medium gain the damping de- 
creases and the phugoid becomes oscillatory and divergent, and for very 
high gain the phugoid becomes aperiodic and divergent. 
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Unlike angle of attack, the velocity feed-back has the principal 
effect on the phugoid. For small gain, the frequency of the phugoid de- 
creases, damping remains about constant, and the phugoid tends to become 
aperiodic, convergent. With another small increase in gain, the phugoid 

becomes aperiodic, divergent; further increase in gain results in greater 

<*• 

divergence. In contrast, the damping and frequency of the short period 
mode are in general only somewhat increased. 

The pitch angle feed-back, the only feed-back designed for stabili- 
zing the motions of the modified Navion, also has a strong effect on the 
phugoid. At low gain, the frequency of the phugoid decreases and the 
damping increases until the motion becomes aperiodic, convergent. The 
short period is effected also, since the frequency increases and the 
damping decreases; this is again a small effect, however, as in the case 
of velocity feed-back. 

Although limited time did not permit an investigation of the effects 
of the feed-backs taken simultaneously, at the magnitudes designed for the 
Ato-2 the separate effects of the feed-backs may be qualitatively summa- 
rized as follows: velocity and pitch angle feed-backs decreased the fre- 

quency of the phugoid of the Navion, thereby changing the frequency to 
that of the A4D-2 phugoid. The angle of attack feed-back was used to 
decrease the frequency and pitch -rate was used to decrease the damping 
of the short period. Finally, both pitch angle and pitch -rate were used 
to increase the damping of the phugoid. 

FLIGHT TEST RESULTS 

As noted in the section on EQUIIMENT AM) PROCEDURE, the in-flight 
simulation of the Al+D-2 was not attempted; some of the reasons therefor 
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are outlined in the preceding section. Flight testing that was done con- 
sisted of individual variations of the feed-back gains, also as mentioned 
previously. Although the aircraft, the variable -stability installation, 
the telemetering station, and associated equipment, were thought to be in 
good working order at the time of the test flight, numerous difficulties 
led to completely unsatisfactory flight test data. 

Problems that arose with the Navion may be enumerated as follows : 
during flight testing of the phugoid mode, it was noticed that engine 
speed varied approximately 150 r.p.m. between high and low speeds of the 
phugoid. It was felt at that time that the added moment from change of 
thrust was causing an alteration of the dynamic characteristics of the 
phugoid. This was particularly noticeable in the testing of the effects 
of the feed-back u; for a small elevator impulse, with resultant small 
amplitudes of oscillations of attitude and airspeed, and with resultant 
small changes in engine r.p.m., the phugoid was a stable oscillatory 
motion. At the same feed-back gain setting, however, and with a larger 
elevator impulse, the phugoid became an unstable oscillation, during which 
time the engine r.p.m. changed considerably. 

Another problem encountered with the Navion was the difficulty in- 
volved in trimming the aircraft prior to each test run. It had been found 
on a previous test flight that the auto-pilot trimming circuit was entirely 
too sensitive and .could not be used to accurately trim the aircraft longi- 
tudinally. Since too little time then remained for changes of trim 
sensitivity, it was decided that the aircraft would have to be re -trimmed 

R 

with the ship's elevator trimming system between each run of the test 
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flight with the auto-pilot disengaged. This was done during the test 
flight, but this method did not provide the cure: the Navion trimming 

system was also found to be very sensitive, and, in addition, manifold 
pressure required to fly the aircraft at the test speed in a level con- 
figuration varied by as much as four inches of mercury. The source of 

j 

this additional trouble was not found. This difficulty in assuming the 
trimmed configuration naturally affected the phugoid test results. 

As will be discussed later, the sensitivity of the pitch channel 
servo amplifier had to be reduced to prevent oscillation or cycling of 
the elevator at the flight test gain potentiometer settings (feed-back 
gains). This reduced the resolution of the elevator servo. The low reso- 
lution of the elevator servo was evident during the tests of the 9 feed- 
back; at low amplitudes of oscillation, there appeared to be little 
response from the elevator servo. During larger amplitudes of oscil- 
lation, the elevator was being moved in what was felt by the test pilot 
to be rather large steps, as indicated by somewhat rough control of the 
aircraft by the auto-pilot and seemingly large step-type changes in 
attitude . 

Additionally, some turbulence was encountered at the test altitude, 
and attempts were made to avoid flying over freshly plowed .fields, and 
the like. Also, communications proved to be unsatisfactory, despite the 
use of a communications relay station between the Navion and the telemetering 
station. Finally, considerable "drop-out" of the telemetered signal was 
encountered at the telemetering station. 
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The result of the above was unsatisfactory flight test data, and 
there was no time left for repeated testing. 

THE VARIABLE-STABILITY INSTALLATION 

The variable -stability installation was built around a three -channel 
auto-pilot, which in turn consists essentially of three servo amplifiers 
and three servo motors used to position the control surfaces. The electric 
stick is used to generate voltage signals, these signals are summed with 
voltage signals from aircraft attitude and rate sensors, and the resulting 
voltage causes the servo motors to position the control surfaces. The 
installation is pictured schematically in Fig. 1. Considering the longi- 
tudinal channel, the servo motor positions the elevator through a bridle 
fastened to the ship's elevator cables, as may be seen in the schematic 
drawing. Fig. 2. Also, since the various sensors and the electric stick 
are excited by a.c. transformers, the signals generated have both magni- 
tude and sense. 

* 

The gain potentiometers provided are used to divide the sensor signal 
voltages, and the portions of the sensor signal voltages summed and fed to 
the amplifier are dictated by the gain pot settings. The gain pot settings, 
in turn, are a function of required gain constants . The algebraic signs 
of the gain constants are determined by the polarity of the exciting trans- 
former, e.g., for the A4D-2 simulation, polarity such that down elevator 
results from an increase in u. 

In addition, a trim circuit was installed. The trim circuit pro- 
vides a voltage of magnitude equal to but of phase opposite to summed 
voltages from the sensors, thereby permitting auto-pilot engagement. 
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To illustrate the elevator-movement sequence, assume that the initial 
input is a voltage generated by stick deflection. The servo and bridle 
move the elevator until the servo follow-up circuit generates a signal 
equal in magnitude and opposite in phase to the stick input signal. As 
the aircraft responds to the elevator deflection, voltages are generated 

Jt 

by the pitch gyro, pitch-rate gyro, airspeed transducer, and angle of 
attack vane. These voltages, after being summed in the pitch channel sum- 
ming circuit, as may be seen in Fig. 3 , are fed through the servo amplifier 
to the servo motor, which deflects the elevator. The elevator augments or 
damps the aircraft response, depending on the phase of the summed a.c. 
voltage fed to the servo motor. 

A detailed description of the variable -stability installation is con- 
tained in Appendix C, including pictures of the installation, descriptions 
of the structural and electrical components, and graphs of the control cable 
forces which were the basis of the structural design. 

Safety features of the auto-pilot consist of the following devices: 
an auto-pilot disengage switch is provided on the left console of the cock- 
pit for the safety-pilot, and a disengage button is provided on the electric 
stick for the test pilot. In addition, an electrical cut-off device is 
provided. If a large error signal is received at the point of input to 
the servo amplifiers, an electrical device breaks the circuit to the 
servos, thereby disengaging the auto-pilot. When this occurs, another de- 
vice disengages completely the gearing between the servo motor and the 
servo pulley; the pulley then rotates freely on bearings only. Finally, 
as noted previously, servo clutches may be caused to slip, i.e., the 
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safety -pilot may manually override the servo motors with the ship's 
controls . 

Structurally, the auto-pilot variable -stability installation was 
found to be satisfactory in all respects except for cable slop in the 
aileron system. Some bending of the aileron servo cable pulley brackets 
was discovered during the static tests; this bending resulted in aileron 
cable slack at very high control force conditions . These pulleys and 
pulley brackets were the original ones used in the Navion aileron follow- 
up system, and because of position, the brackets were not readily ac- 
cessible for structural stiffening. 

AUTO -PILOT LIMITATIONS 

As mentioned previously, the gain constants required for the simu- 
lation were found to be of such magnitude that, for example, the auto-pilot 

.0 o 

had to be designed to furnish about .6 of elevator deflection for a 15 
change in angle of attack. Considering excitation voltage on the servo 
motor (corresponding to about 4 of elevator travel per volt), for a 15 
change in angle of attack, the resulting voltage (from the angle of attack 
vane) which was fed into the summing circuit was about .12 volts. The 
other gain constants, in general, resulted in similar magnitudes of volt- 
ages, with the exception of the pitch angle gain constant, which was an 
order of magnitude smaller. In order to provide such magnitudes of 
voltages, with the available sensor excitation voltages ( + 7*5 to + 30 
volts), and in order to provide gain potentiometer settings that were 
reasonable (towards the middle of the pot indicator scales), 15,000 to 
500,000 ohm end-resistors had to be installed in the individual sensor 



28 . 



I 



circuits in series with the gain potentiometers. (For comparison, gain 
potentiometers and sensor potentiometers have resistances of approxi- 
mately 5000 ohms.) 

In a chain type reaction, considerable quadrature was produced in 
the summing circuit, the trimming circuit could not "interpret" the phase 
shift of the voltage, with the result that the voltage in the summing 
circuit could not be zeroed, and the auto-pilot could not be engaged. 
Quadrature was forthwith adjusted, the auto-pilot was engaged, but with 
any change of gain pot setting or sensor position, the electric cut-out 
disengaged the auto-pilot. Quadrature, then, was found to be a function 
of not only pot settings, but sensor positions as well. In addition, as 
indicated, the electric cut-out could not "distinguish" quadrature from 
a large error signal. 

To relieve the above situation, the emphasis had to be shifted to the 
second objective, investigation of the effects of variations of single 
feed-backs. The end -resistors were removed, but some quadrature effects 
were still found to be present. It was found that, for example, as the 
angle of the angle of attack vane was changed, a small amount of phase 
shift still took place; the explanation was that with the change of re- 
sistance of the angle of attack pot (associated with movement of the vane) 
came also a change in a.c. impedance, causing phase shift. 

To relieve this same situation, and to thereby permit flight-testing, 
the cut-out circuit had to be adjusted to a much less sensitive level, 
and the sensitivity of the servo amplifier had to be decreased to prevent 
cycling of the servo motor. The reduction of sensitivity of the servo 
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amplifier reduced the resolution of the servo motor. The result was a 
barely acceptable longitudinal auto-pilot, considering the small changes 
of aircraft response which occur, for instance, during the phugoid mode. 

To further complicate the picture, the only suitable pressure (air- 
speed) transducer that could be located during auto -pilot installation 
had a non-linear potentiometer. This could be an additional explanation 
for aircraft -auto -pilot behavior during flight testing. 

As mentioned previously, auto-pilot resolution was an additional 
problem. Initially, resolution of the auto-pilot amounted to approxi- 
mately + .5°, i.e., the auto-pilot could position the elevator within 
one -half degree of the commanded deflection. This resolution was found 
to be a function of the excitation voltage on the servo motor, and the 
magnitude of the signal voltage required to close the servo relfiy circuit. 
To improve resolution, the excitation voltage on the servo motor was in- 
creased to the maximum available in the system, + 30 volts. In addition, 
the peck-size of the relay was decreased with use of the servo -amplifier 
throttling control, or in other words, the amount of time that the relay 
remained closed was decreased. This effectively increased the peck-rate 

of the relay. These two improvements increased the resolution to 
o 

about + 1/4 . In an attempt to further increase resolution, a high- 
frequency sinusoidal voltage was superimposed on the relay, or, the 
"dithering" process. The result was unsatisfactory in that the elevator 
cycled at some sub-harmonic of the dither frequency. Finally, the sensi- 
tivity of the servo amplifier had to be reduced to prevent cycling of 
the servo motor, which was believed to be induced by electrical noise 
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in the system; this reduced the resolution of the auto -pilot, but an 
absolute value of the resolution was not thereafter obtained. However, 
considering the .05° of elevator required by the 0 feed-back in par- 
ticular, it may be seen that even the improved resolution is an order of 
magnitude less than that required. 

In summary, the resolution and the quadrature problems were additional 
Indications that a small, auxiliary surface should be used as the moment - 
changer, instead of the Navion elevator. In such a case, voltages would 
have been much larger, and quadrature problems of such magnitude could 
have been avoided. Another approach would have been to design the auto- 
pilot variable -stability system for d.c. summing circuit operation rather 
than a . c . 

CALIBRATIONS 

Quadrature also affected calibrations. With the end -resistors in 
place in the sensor circuits, for accuracy, quadrature had to be constantly 
adjusted during calibration for each range of pot settings. In addition, 
the sensor output voltages, as measured in the summing circuit, were so 
small that "they" could not be measured accurately with the existing volt- 
meters; this necessitated the removal of the end-resistors, calibration 
with other much smaller end -resistors, and a proportionate scaling of 
the calibrated quantities. Final calibrations may be seen in Appendix C; 
they are approximate for no end-resistors in place in the sensor circuits. 
WEIGHT AND BALANCE LIMITATIONS 

According to the Navion Operator's Handbook, maximum permissable weight 
is 2750 lbs., and maximum C.G. position is 30*5$ of M.A.C., gear down. With 
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loaded configuration, including two pilots, full fuel load, and no passen- 
gers, the variable-stability Navion weighs 2750 lbs., and with gear down, 
the C.G. is 29.5# of M.A.C. 

ADDITIONAL DISCUSSION 



In summary, it has been shown in the preceding sections and in the 
Appendices, that with the Navion variable -stability installation, only the 
Navion moment equation was altered, not the lift or the drag equations . 
Further, provisions were made only for modifying two stability derivatives 
of the moment equation, and and for the addition of two pseudo- 

derivatives to the moment equation, and C^. It has been seen that 
with such changes of the Navion moment equation, the characteristic 
equation of the A4D-2 could be simulated, as was verified by use of the 
analog computer. Such simulation meant that the longitudinal dynamic 
characteristics damping and frequency of both modes of motion could be 
simulated. As was shown later, however, steady-state values and ampli- 
tudes of oscillations of the airplane responses could not be simulated. 
This last, coupled with auto-pilot quadrature and resolution problems, 
led to the abandonment of the simulation. 

Although the Navion variable -stability installation has a number 
of problems in the longitudinal channel, this does not detract from the 
value of variable -stability installations in general. The literature 
contains the results of many successful installations, and the real 
value of such an installation is the fact that artificial control per- 
mits changes in the stability characteristics of the test aircraft while 
the aircraft is airborne. Without artificial control, it is impossible 
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to modify existing stability derivatives of a given airplane without 
major and costly structural changes. At the same time., without arti- 
ficial control, it is impossible to add new and additional stability 
derivatives, other than those normally contained in the three equations 
of motion, since there is no choice of what the airplane can "sense." 

With the variable -stability auto-pilot, the outputs of many types 
of sensing devices can be used, and there are practically endless possi- 
bilities in this direction, according to the literature, including Ref. 15 . 
For example, in Ref. k, the investigators used c* and <x feed-backs to modify 
only the characteristics of the short period mode of the test aircraft 
these response feed-backs were used in order not to appreciably change 
the phugoid. Also, u and u feed-backs have been used to only modify 
the phugoid response. As has been discussed in the present report, o<. and 
u feed-backs actually do just what is indicated above, and a quick ex- 
amination of the root locii plots will show that c< and u also have the 
effects described above. 

Another approach to the variable-stability installation theory might 
be changes of the other two equations of motion, lift and drag, although 
this is not believed to be as practical as the moment -changer version. 

Also as indicated in the literature (see Ref. 4), the small, aux- 
iliary surface as a moment-changer is a practical appendage of the variable- 
stability installation. This has been verified by the results indicated 
in the present report. Instead of using a very powerful tool such as the 
Navion elevator to modify aircraft responses, it has been suggested, for 
conjecture, that the elevator trim tab could be used as the moment -changer. 
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Unfortunately, in the case of the Navion installation, this would mean 
additional weight for an additional servo motor, and other auto -pilot 
components, and it has been seen that maximum weight and C.G. conditions 
presently exist. 
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CONCLUSIONS 



In view of the previous, it is concluded that: 

1) Simulation of the longitudinal flight motions of other 
aircraft cannot be successfully accomplished with the 
present Navion variable-stability installation. 

2) Considering simulation in general, unless all three 
equations of motion can be modified, steady-state values 
of aircraft responses and amplitudes of oscillations of 
responses cannot be exactly simulated; however, character- 
istic equations and thus dynamic response characteristics 
such as frequency and damping of the aircraft modes can 

be simulated. 

a) Considering gains of the airplane transfer functions, 
and auto-pilot resolution, a small, auxiliary surface 
should be used as the moment -changer in such type of 
variable-stability installation. Also, the test air- 
craft should have room for weight growth. 

b) Considering gains and numerator terms of the airplane 
transfer functions, in general, the simulator aircraft 
should have approximately the same time constant, 7" > 
as does the aircraft which is to be simulated. 



35 - 



3) In the present configuration the Navion is not suitable for 
in-flight testing of the effects of airplane longitudinal 
response feed-backs. 

4) Concerning feed-back quantities, ot has a primary effect on 
the short period mode, u and -0- have primary effects on the 
phugoid, and -0- has effect on both the short period and the 
phugoid modes . 
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RECOMMENDATION 



Considering present Navion weight, present C.G. location, 
longitudinal auto-pilot limitations, and engine (r.p.m.) 
characteristics, it is recommended that future investigations 
with the Navion variable-stability installation be limited to 
lateral-directional flight testing. 
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NAVION SPECIFICATIONS 


Model * 


North American Navion 


Length 


27.25 ft. 


Wing Span 


33o6 ft. 


Wing Area 


184.2 ft. 2 


Gross Wt. 


2750 lbs. 


C.G. Position 


29.5$ M.A.C. (with Auto -pilot) 


Power Plant 


Continental, six cylinder, 
horizontal-opposed, 205 HP 


Propeller 


Hartzell, variable pitch, . 
constant speed 


landing Gear 


Tricycle, retractable 
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PITCH CHANNEL AUTO-PILOT 
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FIGURE 2 

CONTROL SYSTEM SCHEMATIC 
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PITCH CHAJTT L SU'CIBIG CIRCUIT 
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Terminal number, telemeter a.c. voltage 



APPENDIX A 



THEORY AND ANALYSIS 

MODIFIED NAVION DERIVATIVES AND EQUATIONS 

Given conditions: Clean configuration, power on (1850 rpm) 

h = 6500 ft., 0 = .001957 slugs/ft. 3 
= 120 mph = 176 ft. /sec. 

= .384 (Neutral Point) 

= .095/cleg. = 5»45/rad. 

= - .025/deg. = -1.435/rad. 

= - .15/rad. ( ly) jig = 3055 slug ft. 2 
= - .07/rad. 



From Ref. 11 



Assumption 



Calculations 
C m^ : 



W/S 



C D. 



V, 



N r 



J cK, 



% 



Cm, 



de 



That the derivatives above, which are 'for the 
Ryan Navion, may be applied also to the North 
American Navion; both airplanes are standard 
"blueprint copies" made by two different air- 
craft companies . * 



m c* 

W/S 

% 

CD * 



■ C L * (Xcg - N q ) = 5-M.295 - -384) 
= - .485/rad. 

= 2750/184.2 = 14.93 lbs. /ft. 2 

W/S 14.93 - , ^ 

<1 " l/2( . 001957) (l78) 2 " * 493 

2 C Lc( C Ro( 2( .493) 5.^5 
= = "tt(6.64) .85 = • 333/rad. 



Calculations 



r 

M- 



I CG 



V 



Cl ( k93y 

c D n = c D f + ‘tFaI * 02 5 + V(6.oU).fl5' = - 0401 



r 

A 



m 



2750 



SV 32.17( . 001957 ) ( 176 ) ( 164 . 2 ) 



= 1.35 



l CG 



pSc “ .001957(154.2)5*7 1+1,6 

= lY Jlg - ^ 

= 3055 - 85 .5(1. 39) 2 = 2890 slug ft. 2 



K, 



7 



2 -? -^- 33 . 8«. 8 



l : 


h 


2 *y* 

/J- 


• 

“u 


° m u 


i 0 


D eff 


C;D eff 


± n 

D o 


1 

L eff 


Cr 

L eff 


= C L 
L o 



55^5 

2 . 33*8 , 

5TT TJTrT 



by - h r 



Characteristic equation. Modified Navion: Dividing all coefficients 

4 



c^ A ^ + c 3 X -^ + c 2 K 2 + Ci X +c Q =o 



where 



C4 = 1 



= fc 

i> r L 



C m Go + C m - h( Z° L- + C 
m d0 2 m d<X ' 2 



Deff^] 



= .15 % - .07 - .05(2.725 + .o4oi)j 

= 5*30 + 15*75 Kg 



A 2 



sec . 



C 2 - [% ff ( Si* * c -ae G2 ' *1 £WCd k ' ° l) 

+ ^^ C md0 G2 ^ + Cm^ G 1 + K 4 C m&JJ 

- ■■ 65 bfr.(fe 27 [• <A01 { - - 07 • - 15G2 - •°5(2.T25)| 

. .05(.493) (-333 - .493) + 2 -7 2 5( " ^2) " - 48 5 & 1 
2 

- 1.435%J 

= 9.88 + 15.74% + 32.27% + 15.73% 

C 1 = 3^.3 [ C D eff (^J* C mde G 2 + C m ^ G x ) + ^L^L eff (C md0 G 2 



+ C m de ^2) - ^ ^ + < JCl 



= -1 
.05(2.46) 



.493 ) 

2 



2 ( 



-.150, 



L - W-H 

0401^2. 72 5( -.15)02 " ‘485G^ + (_. 

- .07 - (- .15)G 2 J - ^ (.333)03- i.435K4(. 0401+2. 725) 

= .4o8 + .468% + 1.284% - 1.942% + 32.25% 

C = i — 4 1 C L Ct C m G-i - C L C m Go] + KkC m Ct, __ Ell* 

-hr \ 2 ~ Leff m <* 1 g-‘ * “u 3 4 ra ^LDeff 2 

- - c L )V 



- ' 0m 32 [-*93(. | «)6 1 • 2 - 725 Vs] 

+ ^ A3 ^ K4 [(.0401)2.725 - • 4 |3(.333 _ .493)] 

(.05)3.32 L 2 J 



= .355 + 1.051 % - 5.81 % + 1.285 % 
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A4D-2 derivatives and equations 



From Ref. 12, the A4D-2 derivatives are as follows 
Given conditions: Clean configuration, power on 

h 



From Ref . 12 : 



= 6500 ft. V Q = Mach .2 = 218 ft. /sec. 



C T = 3.62/rad. 
L <* 






m 



c< 



= - .3265/rad. 
= - .ll+5/rad. 



= - 1.090/rad. = 

ra D* 



in 






m 



3 .263/rad. 



^ lv)' 

b l ~'m 



b(ac) 

2V 



“u 

% 

G 

CG 



= 0 

= 17,600 slug ft. 2 , S = 260 ft. 2 
= 10.8 ft., W = 10,000 lbs. 

= 25$ MAC, b = 27.5 ft. 



Computation: 



"Dr 



"D 



<* 



■'min 



_ W/S _ 10 _ OpO 

q. l/2( .001957) (218) 2 260 '° 2 ° 

i .059, e ± .572, A = 2.91 



± C T 



■mn 



+ (cL ) 2 

-rrA e 

urn 



- .059 + ^.*2.91( .572) 1 * 190 

C D = 2Cl,CLpt = 2(;828)3-62 = i.iLy/rad. 

D « jta p jr(2.9i).572 



r 



r 



m 



10, 000 /g 



pSV .001957 x 260 x 218 



=2.81 sec. 



Computation: 



C m 

m d<* • 


Cm dc< 


= r. c 1.090(10.8) = . 

m D*'2 vr 2(218)2.81 


.00961/rad. 


Cm de : 


J 

Cm d0 


£ -3.263(10.8) 

- • 2VA 2(218)2.81 


.0288/rad. 


ro 

•> 


h 2 


- *v - 17 ,600 ^ £,4. 2 

- l 0 ;ooo/ g -56-6 ft. 




A : 


A 


m 10,000/g £ 

~ ~psc ~ .001957(260)10.8 " 50 




h : 


h 


2 AC 2 (56.6)„ 

if = 5^tf nSSj 2 - - 01715 



Characteristic equation, A4D-2: Dividing all coefficients by -hf 

' , 4 t,3 i , 2 t . ' 

Ci| }\ + C^ + Cg A + A + C Q = 0, where 

Ci| = 1 



c 3 - hT- [■ h(C » + -4^ ) + Co a« * C “. 



de 



= (. 01715 )C2.bl) 



m dw + Cm d(9^ 



+ ■ u ■ C 



[- . 01715 ( .19 + 1.81) - .00961 - .0288 
h [ h ( \ - c o - + c d (c 

+ C™ I 
m <x 

= . 01715(2. 81 ) 2 [-01715 ^.414(1.147 -.828) -.19(1.81)1 

+ .19(- .00961 -.0288) + l.8l(- .0288) - . 3.45^1 
= 1.536 



= 1.508 



A 5 



I 



C- - il_ 

°i “ h r 



■3 [ c d(^2. *c mc< ) + EHc^ i^- c L | ) 



= ( .01715) (2.61)5 
+ .0288( .319)|| 

= .0968 

2 , 



[•19(1. 



8ix-.0288 -.1U5) + .4 i 4\(.828)(-. 00961) 












GAIN CONSTANTS 

Equating the constants of the two characteristic equations : 
( C 0a 4D-2 = ( C 4VoD.NAV. = 1 

% = 1.508 = 5.30 + 15.75% 

c' 2 = 1.536 = 9-88 + 15-74% + 32 . 27% + 15-73% 

% = .0968 = .408 + .468% + 1.284% - 1.942% + 32.25% 

% = .0464 = .355 + 1.051% . 5.81% + 1.285% 

Solutions of these equations yielded: 

% = - .0407, dimensionless 

% = - .240, sec. 

% = .0465, rad. 

% = .0033, dimensionless 
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APPENDIX B 



ANALYTICAL STUDIES 

COMPARISON OF DYNAMIC CHARACTERISTICS 

Factorization of the characteristic equations yielded the following 
AUD-2, and MOD. NAV.: 

Rrngoid X 2 + .035lX + .0325 = 0 

X phu = - .01755 + .I8i 
T l/2 = 39 - 5 sec. P= 34.9 sec. 

Short Period X 2 + 1.473 X + 1.452 = 0 
X sp = - .737 + .9541 
T1/2 = -94 sec. P= 6.59 sec. 



NAVION : 

Hiugoid 



X 2 + .0212 X + .0366 = 0 

T , = 62.6 sec. P= 32.9 sec. 

i l/2 



Short Period X + 5-28 X+ 9-72 = 0 

T1/2 = *26 sec. P= 3*78 sec. 

COMPARISON OF STEADY-STATE VALUES 

The steady-state denominator term common to each aircraft transfer 
function : 





Cl) eff 


1/2 (C D - C L ) 


C L/2 


A,= 


C Leff 


C 4<x 

2 


0 




0 


c °v 


0 




c L 2 c m 








2 


B1 





Steady-state numerator terms: 



(»u) 



S .S . 



0 



- S, 1 * 



1/2(C e - C L ) 



2 

C nu 



c L c L 

m S e e I r L * 



("*> 



s .s , 






eff 



J eff 



0 

0 

S« ^ 



C L 

2“ Cm s/ e 



0 



0 



C L/2 



0 



0 





Cr) eff 


1 / 2 (Cp^ - c L )' 


0 


( ) S . S . = 


C Leff 


Ct 

h 

2 


0 




0 


Cm o«. 


~ c m 




= " Cm s a ^ 


r* - & 

2 2 * 


- c L ) 



A4D-2: 

Evaluation of steady -state values for .05 rad. 

(down eldvator) step forcing function yielded: 
u s.s. = + *2^6 ^Computer: + .255 J 

0<s.s. = - -H3 rad. = - 6.45° [^Computer: - .114 radQ 

9 S . S . = - «069 rad. = - 3*95° [^Computer: - .070 rad/j 
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MODIFIED NAVION: 



For the 


same forcing function as 


above (including feed-backs) 


u s .s. 


= 6.15 


^ Computer : 


5.8 ] 


^s .s . 


= - 1.1 rad. 


[^Computer: 


- 1.0 rad/J 


®s .s . 


= - 1.35 rad. 


^Computer : 


- 1.3 radfj 



COMPUTER (DYNAMIC STUDY) 

A 4 D- 2 : Computer equations: 

C Qi 

Lift : u+ ~jr u +_ 2J“(^D 0 ^ ” + = ^ 

C L Ct 

Drag : c*+ -p u + ^7 °< - © =0 

Moment : 9 + (’^ )«+ ( 1 ^ )* + &£)© + (~° m 6 e )L = 0 

h r*- hr hr' hjcir 



Damping and Period: 

Riugoid Tj/ 2 = ^ 0 *° sec *> P =35*0 sec. 

Short Period t i/2 ® sec., P i 8.0 sec. 

Transient values: For .05 rad. <S e (up elevator) 

step input, amplitudes of initial 
oscillations are as follows: 
u = - A5 



°*sp - 5*5 

°^phu = 7 ' 8 ° 

© = 19. 5° 
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Transient normal acceleration: 



MOD. 

Lift 

Drag 

Moment 



W(n - 1) = g a 

/ X v 2 • V 

or (n - l)g = a = — , and © = ^ 

then © = g^ n ^ 



vd 



or n = — + 1 . 



From the computer, 

© = .08 rad. /sec. initial amplitude of oscillation 



V© 

n = g 



+ 1 = 



during the short period; 

219(.08) 

32.2 



+ 1 = 1.5 g's 



NAV . : Computer equations : 

u + .0297u - .0592 <* + .18259 = 0 

ck + . 3655 u + 2.02 o< -9 = 0 

© + .730u + 5.32<x - .641 tx, + 1.038 oc + .0518© + 2.22© 

- 3.78© - 15.75 Ks & s = 0 
Damping and Period: 

Riugoid ^l/2 = 39*0 sec., P = 35*5 sec . 

Short Period T^y^ = *5 sec., P = 4.8 sec. 

Transient values: For .05 rad. & e (up elevator) step input , 

amplitudes of initial oscillations are 
as follows: 
u = -10 
c< = +83° 
n = 6.8 g's 
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For .05 rad. 6 e (up elevator) Impulse 



forcing function, initial oscillations: 

= 12° [Function of feed-backs^ 
u = - -75 
«sp = 16 ° 

? o 

oi . =4 

phu 

n = 3.2 - 4.0 g's 
COMPARISON OF TRANSFER FUNCTIONS 

The denominator term common to all of the airplane transfer functions 

(not including auto-pilot dynamics) may be written as: 

4 

A (S) = - hi" (S - A x )(s - b 2 ){ S - * 3 )(S - * 4 ), Where 

^2 ^ are the roots of the characteristic equation 

(A4D-2 or MOD.NAV.). Numerator terms may be written as follows: 



N u - + c m 



V ° £ 



4, r 



[ s + C L* C L 1 

2 L ^ \ 



N * - 






[■ 



C mc 4 e r 2 [s 2 + ^S + l/2(^) 2 ] 



w o - - S 



e S t r 2 [s 2 + |( C D + S +23^ (C D C V C L CP 



+ C L ) 



Now 



U / X 

examine a typical transfer function, say 7\S): 



U_ (s) = . h . P 

- 2hr" 



r 

L +2TC Vo 



[_ (s -A x ) (s -A 4 )] 
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Evaluation of this transfer function gives: 

1.941 Ls + 2.99] 

MOD. NAV-: = 



[ (s - Xj) . . . . (s -^ )] 



.492 [s + .465] 

a 4 D -2 : = JJs - A 1 )....(S 







ibis may also be written in servo-mechanism notation as follows: 



J*( S ) „ K[s 2 , g ^nS ^n 2 3 „ here 
[(s - * x ) (s -a 4 )] , 

K is the gain, 5 is the damping ratio, and u) n is the natural frequency 
of this minimum phase term (or second order lead term). 



Evaluation of the transfer function gives : 



MOD. NAV.: 
A4D-2 : 



- 15-75 [s 2 + 2 (. 0576 )(-g 58 )S + (. 258 ) 2 ] 

= [(s -Ai> (s - \) j 

- 2.41 [s 2 + 2(.l62)(.208)s + (.208) 2 ] 
[(S - *i) (S -* 4 ) J 



And examine 




9 , 

Tl < s > : 

[s 2 + i (Cj^+ Jw)s + (^1. - C L C Ew + C l 2 )J 

[ (S- * x ) (5-^4) ] 



b6 



- 15-75 Ls 2 + 2(3-58)(.286)S + (.286) 2 ] 
MOD. NAV.: = [( S - (S - * 4 ) ] 

-2.41 [s 2 + 2(2.17 )(.i 64 )s + ( .164 ) 2 J 

A4D-2 : = £( S _ A x ) (S - *4)] 



ROOT LOCUS STUDY OF FEED-BACKS 

In a study of the feed-back quantities used to simulate the A4D-2, 
the characteristic determinant for the Navion must first be examined. 
With angles in radians, and in real time, the determinant is: 

o< 0 

l/2(c^ - c L ) p 

J*p> - r P 

+ Cma<* r p) ^ Cm d0 “ 11 ^ P^P 

With the addition of a u feed-back in the moment equation, the 
moment equation takes the following form: 

•S V u + (c -v * + »r p ) rp8 = - c m& s e 

Then the characteristic equation takes the following form: 

- hr 4 u- Aq)( ^ - a 2 )U- a 3 )U -H) 



U 

< c d + r P ) 

C L 



S 



: m , r _ 
2 



D 



c< 



C L* 

2C D^ 



= 0, 



where A 4 are the roots of the characteristic equation. 

Such roots naturally are modified by the root which is part of the second 
term of the equation. The equation may be further simplified: 
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a( X - X x ) 



(X -K k ) 



- k 3 = o, 




o 

where a = - h T /C m Cp 



Cp , and is a positive number. Then: 




! d^ navion 




Also, the roots of the characteristic equation of the Navion are: 

4>hu - - -0106 t .1941 
X = - 2.64 + 1.671 

s .p. 

where "phu" indicates phugoid, and "s.p." indicates 
short period. 

Refer to the root locus diagram, Fig. B2. For positive, or 
down elevator for increasing u, examine the 0° condition: it may be 

seen that damping and frequency of the short period mode are somewhat 
increased (period and time to damp to one-half amplitude, P and T \J2> 
are decreased), but the characteristics of the phugoid mode are altered 
appreciably. The phugoid oscillation rapidly becomes unstable, then 
becomes aperiodic and divergent. 

In like manner, feeding back o<, 9, and 9, singly and in order, 
results in the addition of the following terms to the original Navion 
characteristic equation: 
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o< 



K iS/ 






A = - .0147 + .2581 

2 , Q 

a = - hJ“/C m j 0 condition. 



0 : 



+ KgCnj. 



r 2 ^ 



2 1 
1 + r(°D 



+ c H>i 



CdCl 

2T^ 



<x 




C l) 



6 



A = 0, - 2.01, + .04 



■» = - hJ? / C m s J 

'V.l/l '' 2 



Q°condition 
1 Ct 

+ jC (Cj)+ 2 — )A 



C D C L« 
— vg" 

2 J" 



c L 

- (C Dp< - C L ) 



A = - 2.01, + .04; 180 ° condition 
For the results of the root locus studies, refer to Figs. B2 
through 5* For a summary of results of Appendix B, refer to Tables 
I and II. 
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liary Circuit for Elevator Deflection 
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APPENDIX C 



VARIABLE-STABILITY INSTALLATION 

The selection of the Minneapolis -Honeywell, USAF TYPE E-12, Auto- 
pilot imposed some limitations on the force and speed of response avail- 
able for control surface movement. Flight test results, presented in 
Figs. Cl and C2, showed that the force requirements of the vehicle for the 
variable-stability installation, Princeton University's North American 
Navion, (Fig. C3)> were well within the l60 in. lb. shaft-torque output of 
the servo motors. The speed and frequency response characteristics were 
studied by another group working on the installation and are discussed in 
Ref. 9. 

Blueprints of an auto-pilot installation designed for the Navion were 
obtained from Lear, Inc. From this basic design two brackets were made 
of l/8", 24-ST, aluminum sheet and bolted to the partial bulkhead at the 
after end of the baggage compartment. Across these vertical brackets 
a l/8", 24-ST, aluminum plate was bolted such that the plate lay in a plane 
approximately 90° to the ship's control cables. To support the servo motors 
a four-sided box, about eleven inches square, was fabricated of l/8", 24 -ST, 
aluminum sheet and 3/4" by l/8", 17 -ST, aluminum angle extrusions. To each 
of three sides of the box a servo was bolted, the bolts extending through 
aluminum angle stiffeners on the inner faces of the box sides. The remain- 
ing side of the box was bolted to the plate and bracket combination. Ad- 
ditional stiffening was supplied by l/l6", 24-ST, aluminum sheeting across 
the top and bottom of the box and two small l/8", 24-ST, sheet metal fillets 



Cl 



bolted to the vertical brackets and the partial bulkhead as lateral supports. 
The entire servo mount was over-designed to preclude structural vibrations; 
the completed unit is shown in Fig. Cb. 

With the servos mounted vertically the capstans, or motor-driven 
pulleys, werg nearly parallel with the ship's rudder and elevator cables. 

The bridle cable from the elevator servo, made of 3/32" aircraft control 
cable, was fastened to the capstan and fed directly aft about five feet to 
join the elevator cables at a small angle. The l/8" rudder bridle cable 
was wrapped around the capstan and fed directly aft on the right side; a 
two in. fiber pulley was mounted on the under side of the sheet metal box 
to aline the bridle cable with the left rudder cable. The aileron bridle, 
of 1/8" cable, also was wrapped around its capstan, and the free ends of 
the bridle were fed forward through existing cable access holes in the 
bulkhead; a two in. fiber pulley was necessary for alinement on the right 
side. The 1/8" bridle cable was then mated to the 3/32" cross-control 
cable, the latter having been disconnected from the rudder cables. Standard 
aircraft cable thimbles and turnbuckles were used in all cable connections . 
The rudder and elevator bridles were attached to the ship's cables using 
sheet metal c -clamps. The bridle cable was united with the clamp by 
a 3/l6" bolt passing through a thimble and the clamp; the clamp was tight- 
ened around the control cable with three more bolts. After completion of 
the bridle and cable matings, turnbuckles were adjusted for 25 # tension 
in all bridle cables and 3 0# tension in the ship's cables. Servo clutch 
slip levels were adjusted to allow overpowering of the auto-pilot with 
reasonable control forces. The electrical servo stops were set to prevent 
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full deflection of any control surface in response to auto-pilot signals. 

A schematic diagram of the bridle -control cable arrangement is shown 
in Fig. 2. 

The conventional yoke and rudder pedals were removed from the right 
side of the cjockpit. Two spring-loaded rudder pedals were mounted on the 
floorboards; a small cable attached to the pedals was used to drive 
a 36 O 0 potentiometer. The electric stick, a spring -centered flight con- 
troller with adjustable fluid damping, was bolted to a pedestal for appro- 
priate height. The flight controller had no provision for stick force 
inputs; it had a small break-out force that was not objectionable. In 
addition to the above, Fig. 5 shows the auto -pilot trim box as it was 
mounted adjacent to the flight controller. 

Other modifications to the cockpit included the various auto-pilot 
and telemetry power switches mounted to the left of the instrument panel 
and the feed-back gain potentiometers grouped between the two front seats 
(see Fig. C6). The gain potentiometers are wire-wound potentiometers of 
approximately 5000 ohms each. With these feed-back controls eight air- 
craft responses ( jZS, ty, 6-, ■©, Q ,o< and u ) may be fed into the auto-pilot 
to vary the response, and hence the dynamic stability, of the aircraft. 
Control surface travel can be varied by either the stick and rudder pedal 
gain potentiometers or by the servo follow-up gain controls. 

The major portion of the variable -stability installation was located 
in the baggage compartment just aft of the rear seat. The top shelf of a 
table was devoted to terminal strips that gave accessibility to the wiring 
of the system. These strips made it possible to monitor all the sensor 
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and excitation voltages; also, the sense of control deflection, for example 
down elevator for increasing u, could be changed at the terminal board 
merely by interchanging the transformer leads exciting the signal generator. 
The power to excite the primaries of the four multiple -wound transformers came 
from an invej^ter which converted the aircraft 28 volt d.c. supply to 115 
volt, single phase, 400 cycle a.c. 

The electronic portions of the auto-pilot were placed on a lower shelf 
along with two rate gyros and the pitch -roll attitude gyro; three more rate 
gyros were mounted on the floor of the baggage compartment (Fig. C7)* Al- 
though only three rates of response were measured a separate gyro was re- 
quired for each channel using the rate signal. The directional coupler 
feed-back loop and the control stick and co-ordinated turn capabilities 
of the auto -pilot were eliminated. An electronic safety device was ad- 
justed to disengage the auto-pilot, by releasing the servo clutches 
allowing the capstans to rotate freely, if an error voltage of approxi- 
mately three volts suddenly appeared as an input to any servo amplifier. 

A disengage switch on the control stick was wired also, along with another 
disengage switch, for the safety pilot, adjacent to the main power switch. 

To record in-flight data an ASCOP pulse width frequency modulating 
telemeter station was used. The 15 channel filter conversion unit and 
the rotary sampling switch and transmitter unit were mounted in the 
baggage compartment (Fig. C4). The filter conversion unit allowed sampling 
of fqurteen signals, reserving one filter for reference voltage. The 
rotary sampling switch operated at 20 cps and could sample U3 quantities 
plus two zero readings; hence a telemetry terminal strip was made to 



parallel switch channels so that each of the 14 filter outputs could be 
sampled three times per cycle. Additional data on the internal components 
of the auto-pilot may be found in Ref. 9* 

To accommodate the remaining attitude sensors two booms were mounted, 
one at each ying tip rib as shown in Figs. C8 and C9* The booms, con- 
structed for previous testing, placed the sensing devices about four feet 
in front of the leading edge, ahead of the major pressure disturbances 
around the wing. Dynamic and static pressure lines extend from the 
pressure probe to a differential pressure transducer mounted inside the 
wing. Two Giannini yaw vanes were orientated 90° apart on the left probe 
to serve as sideslip and angle of attack vanes; vane mutual interference 
was considered beneath the resolution of the system. Since two position 
potentiometers were actuated by each vane only one angle sensor was re- 
quired to supply signals to both the yaw and roll channels. 

To complete the installation it was necessary to calibrate the entire 
system, that is to determine sensor signal voltages, for various gain 
potentiometer settings, throughout the operating ranges of the various 
attitude and rate sensors. Also it was necessary to determine control 
surface deflection for given input voltages . Then, to enable reading of 
the telemeter traces, ratios of a.c. sensor signal voltage to d.c. telemeter 
signal voltage had to be established for each of the quantities to be 
recorded. Results of the longitudinal channel calibrations are shown 
in Figs. CIO and Cll. 
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